Critical behavior of n-vector model with quenched 
randomness 

J. Kaupuzs 

Institute of Mathematics and Computer Science, University of Latvia 
Rainja Boulevard 29, LV-1459 Riga, Latvia 
e-mail: kaupuzs@latnet.lv 

Abstract. We consider the Ginzburg-Landau phase transition model with 0(n) 
symmetry (i.e., the n-vector model) which includes a quenched randomness, i.e., a 
random temperature disorder. We have proven rigorously that within the diagrammatic 
perturbation theory the quenched randomness does not change the critical exponents at 
n — > 0, which is in contrast to the conventional point of view based on the perturbative 
renormalization group theory. 
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1. Introduction 

The phase transitions and critical phenomena is one of the most widely investigated 
topic in modern physics. Nevertheless, an eliminated number of exact and rigorous 
results are available, and they refer mainly to the two-dimensional systems |J and 
fractals || . Rigorous results have been obtained also in four dimensions [|3| based on an 
exact renormalization group (RG) technique @]. The RG method, obviously, provides 
exact results at d > 4 (where d is the spatial dimensionality), but this case is somewhat 
trivial in view of critical phenomena. In three dimensions, approximate methods are 
usually used based on perturbation theory. 

Here we present a particular result obtained within the diagrammatic perturbation 
theory. The Ginzburg-Landau phase transition model with 0(n) symmetry (i.e., the 
n-vector model) is considered, which includes a quenched random temperature disorder. 
The usual prediction of the perturbative RG field theory p, §, [7| is that, in the case 
of the spatial dimensionality d < 4 and small enough n (at n — 1 and n — > 0, in 
particular), the critical behavior of the n-component vector model is essentially changed 
by the quenched randomness. Here we challenge this conventional point of view based 
on a mathematical proof. We have proven rigorously that within the diagrammatic 
perturbation theory the critical exponents in the actually considered model cannot be 
changed by the quenched randomness at n — > 0. 
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2. The model 



We consider a model with the Ginzburg-Landau Hamiltonian 



H/T = J [(r + v^/(x)) <f?(x) + c(V<p(x))' 



dx 



+ uV 1 ^(ki^k?) Wki+k 2 ^(k 3 )^(-ki - k 2 - k 3 ) (1) 

i,j,ki,k2,k3 

which includes a random temperature (or random mass) disorder represented by the 
term ^Ju f (x) (p 2 (x) . For convenience, we call this model the random model. In 
Eq. (|I|) <p(x) is an n-component vector (the order parameter field) with components 
<fii{x) = V~ l l 2 J2k ^(^e* 1 ™, depending on the coordinate x, and /(x) = V" 1 ^ 2 J2k /ke 2kx 
is a random variable with the Fourier components /k = V~ l l 2 J /(x)e _ikx (ix. Here V 
is the volume of the system, T is the temperature, and y?«(k) is the Fourier transform 
of ipi(x). An upper limit of the magnitude of wave vector fc is fixed. It means that 
the only allowed configurations of the order parameter field are those corresponding to 
<£>i(k) = at k > k . This is the limiting case m — > oo (m is integer) of the model 
where all configurations of <p(x) are allowed, but Hamiltonian (HD is completed by term 
E(k/k ) 2m \ ^(k) | 2 . 

i,k 

The perturbation expansions of various physical quantities in powers of the coupling 
constant u are of interest. In this case n may be considered as a continuous parameter. 
In particular, the case n — > has a physical meaning describing the statistics of 
polymers 0. 

The system is characterized by the two-point correlation function Gj(k) defined by 
the equation 

(^(k)^(-k)>=5 w G i (k) . (2) 

Because of the 0{n) symmetry of the considered model, we have Gj(k) = G(k), i.e., 
the index i may be omitted. It is supposed that the averaging is performed over the 
/(x) configurations with a Gaussian distribution for the Fourier components /k, i.e., 
the configuration {/k} is taken with the weight function 



P({h}) = Zi 1 exp (-p(k)|/ k |j 



(3) 



where 

Zi = /exp(|-X:6(k)|/k| 2 J D(f), (4) 

and b(k) is a positively defined function of k. Eq. ([!]) defines the simplest random model 
considered in (according to the universality hypothesis, the factor \Ju does not make 
an important difference). Our random model describes a quenched randomness since the 
distribution over the configurations {/k} of the random variable is given (by Eqs. (^) and 
(U)) and does not depend neither on temperature nor the configuration {^(k)} of the 
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order parameter field. More precisely, the common distribution over the configurations 
{/ k } and {<#(k)} is given by 

P({/ k }, to(k)}) = P({/ k }) x Z;\{U}) exp(-H/T) , (5) 

where Z 2 ({f\ l }) = J exp(—H/T)D((p) and H is defined by Eq. ([I]). For comparison, 
the common distribution is the equilibrium Gibbs distribution in a case of annealed 
randomness sometimes considered in literature. 



3. A basic theorem 



We have proven the following theorem. 

Theorem. In the limit n — > 0, the perturbation expansion of the correlation 
function G(k) in u power series for the random model with the Hamiltonian ([!]) 
is identical to the perturbation expansion for the corresponding model with the 
Hamiltonian 



H/T = J [r ^ 2 (x) + c(V^(x)) s 



<ix 



+ uV 1 y?i(ki)y2;(k 2 ) u kl+k2 V2 j (k 3 )v9 i (-ki - k 2 - k 3 ) (6) 

ij',ki,k2,k3 

where w k = w k - \ (\ / k 

For convenience, we call the model without the term y/u f(x) y? 2 (x) the pure model, 
since this term simulates the effect of random impurities ||. 

Proof of the theorem. According to the rules of the diagram technique, the formal 
expansion for G(k) involves all connected diagrams with two fixed outer solid lines. In 
the case of the pure model, diagrams are constructed of the vertices >«vw*< ; with 
factor — uV^Uk related to any zigzag line with wave vector k. The solid lines are related 
to the correlation function in the Gaussian approximation Go(k) = 1/ (2ro + 2cfc 2 ). 
Summation over the components y^(k) of the vector ^(k) yields factor n corresponding 
to each closed loop of solid lines in the diagrams. According to this, the formal 
perturbation expansion is defined at arbitrary n. In the limit n — > 0, all diagrams 
of G(k) vanish except those which do not contain the closed loops. In such a way, for 
the pure model we obtain the expansion 

G(k) = k ,k + krT^-k + ... . (7) 

In the case of the random model, the diagrams are constructed of the vertices 

>• < and ••■•< . Besides, it is important that only those diagrams give a 

nonzero contribution where each dotted line is coupled to another dotted line. The 
factors uV^ 1 (j / k correspond to the coupled dotted lines and the factors — uV^u^ 
correspond to the dashed lines. Thus, we have 



G(k) 



ki ^k + k_ 



+ ... . (8) 
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In our notation the combinatorial factor corresponding to any specific diagram is not 
given explicitly, but is implied in the diagram itself. In the random model, first the 
correlation function G(k) is calculated at a fixed {/k} according to the distribution 
Z^ 1 ({/k}) exp(— if/T) (which corresponds to diagrams where solid lines are coupled, 
but the dotted lines with factors — y/u V^ 1 ^ 2 /^ are not coupled), performing the 
averaging with the weight (S) over the configurations of the random variable (i.e., the 
coupling of the dotted lines) afterwards. According to this procedure, the diagrams of 
the random model in general (not only at n — > 0) do not contain parts like , 

^•l! > > etc- Such parts would appear after the coupling of dotted lines only 

if unconnected (i.e., consisting of separate parts) diagrams with fixed {/k} would be 
considered. 

Thus, in the considered random model, the term of the I— th order in the 
perturbation expansion of G(k) in u power series is represented by diagrams constructed 

of a number M\ of vertices >■ <C and an even number M 2 of vertices ••••< 

(i.e., M 2 /2 double-vertices >• <C ), such that I = Mi + M 2 /2. In the pure 

model, defined by Eq. ([]), this term is represented by diagrams constructed of / vertices 
>**w*< . It is clear and evident from Eqs. (0) and @ that all diagrams of the random 
model are obtained from those of the pure model if any of the zigzag lines is replaced 
either by a dashed or by a dotted line, performing summation over all such possibilities. 
Note that such a method is valid in the limit n — > 0, but not in general. The problem is 
that, except the case n — > 0, the diagrams of the pure model contain parts like , 

■ ^8c£ ' e ^ c - ^ a ^ ^he depicted here zigzag lines are replaced by the dotted 
lines, then we obtain diagrams which are not allowed in the random model, as it has 
been explained before. At n — > 0, the only problem is to determine the combinatorial 
factors for the diagrams obtained by the above replacements. For a diagram constructed 

of Mi vertices >» <C and M 2 vertices ••••< the combinatorial factor is the number 

of possible different couplings of lines, corresponding to the given topological picture, 
divided by Mi!M 2 !. 

It is suitable to make some systematic grouping of the diagrams of the random 
model. The following consideration is valid not only for the diagrams of the two-point 
correlation function, but also for free energy diagrams (connected diagrams without 
outer lines) and for the diagrams of 2m-point correlation function (i.e., the diagrams 
with 2m fixed outer solid lines, containing no separate parts unconnected to these lines) . 
It is supposed that at n — > the main terms are retained, which means that the free 
energy diagrams contain a single loop of solid lines. We define that all diagrams which 
can be obtained from the z— th diagram (i.e., the diagram of the z— th topology) of the 
pure model, belong to the z-th group. The sum of the diagrams of the z— th group can 
be found by the following algorithm. 

1. First, the z-th diagram of pure model is depicted in an a priori defined way. 

2. Each vertex >^~^< is replaced either by the vertex >> <C , or by the 

double-vertex >» •< , performing the summation over all possibilities. Besides, all 
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vertices >> <C and ••••< and all lines are numbered before coupling, and all the 

distributions of the numbered vertices and lines over the numbered positions (arranged 
according to the given picture defined in step 1 and according to the actually considered 

choice, defining which of the vertices >ww< must be replaced by > <C and 

which of them must be replaced by >• •< ) are counted as different. Each specific 

realization is summed over with the weight 1/(M 1 \M 2 \). 

3. To ensure that each specific diagram is counted with the correct combinatorial 
factor, the result of summation in step 2 is divided by the number of independent 
symmetry transformations (including the identical transformation) Si for the considered 
i-th diagram constructed of vertices >*»w~< ; where the symmetry transformation 
of a diagram is defined as any possible redistribution (such that the outer solid lines are 
fixed) of vertices and coupled lines not changing the given picture. Really, the coupling 
of lines is not changed if any of the symmetry transformations with any of the specific 
diagrams of the i— th group is performed, whereas, according to the algorithm of step 2, 
original and transformed diagrams are counted as different. 

It is convenient to modify step 2 as follows. Choose any one replacement of the 

vertices >^^< by >> <C and >• •< , and perform the summation over all 

such possibilities. For any specific choice we consider only one of the possible M\\M 2 \ 

distributions of the numbered Mi vertices >■ ■< and M 2 vertices ••••< over 

the fixed numbered positions, and make the summation with weight 1 instead of the 
summation over M]\M 2 \ equivalent (i.e., equally contributing) distributions with the 
weight l/(Mi!M 2 !). 

Note that the location of any vertex >■ <C is defined by fixing the position 

of dashed line, the orientation of which is not fixed. According to this, the summation 
over all possible distributions of lines (numbered before coupling) for one fixed location 
of vertices (as consistent with the modified step 2) yields factor g M i4 M 2/ 2 _ The i- 
th diagram of the pure model also can be calculated by such an algorithm. In this 
case the summation over all possible line distributions yields a factor of 8', where 
I = M1 + M2/2 is the total number of vertices >~--~< in the i— th diagram. Obviously, 
the summation of diagrams of the i— th group can be performed with factors 8 l instead of 
8 A/l 4 M2//2 , but in this case factors §wV" -1 (j /k | 2 ^> must be related to the coupled dotted 

lines instead of uV~ l /k | 2 ^>. In this case the summation over all possibilities where 
zigzag lines are replaced by dashed lines with factors —uV~ x u\c and by dotted lines with 
factors \uV~ l /k | 2 ^>, obviously, yields a factor uV^ 1 (— Uk + ^(\ fu | 2 )) = —uV~ x u-k 
corresponding to each zigzag line with wave vector k. Thus, the sum over the diagrams 
of the i-th group is identical to the z-th diagram of the pure model defined by Eq. (|6]). 
By this the theorem has proved. 

4. Remarks and conclusions 

The theorem has been formulated for the two-point correlation function, but the proof, 
in fact, is more general, as regards the relation between diagrams of random and pure 
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models. Thus, the statement of the theorem is true also for free energy and for 2m-point 
correlation function. 

According to the proven theorem and this remark, at n — > the considered pure and 
random models cannot be distinguished within the diagrammatic perturbation theory. 
Thus, if, in principle, critical exponents can be determined from the diagrammatic 
perturbation theory at n — > 0, then, in this limit, the critical exponents for the random 
model are the same as for the pure model. We think that in reality correct values 
of critical exponents can be determined from the diagrammatic perturbation theory, 
therefore the quenched random temperature disorder does not change the universality 
class at n — > 0. This our conclusion correlates with results of some other authors. 
In particular, there is a good evidence that the universality class is not changed by 
the quenched randomness at n = 1. It has been shown by extensive Monte-Carlo 
simulations of two-dimensional dilute Ising models || that the critical exponent of the 
defect magnetization is a continuous function of the defect coupling. By analyzing the 
stability conditions, it has been concluded in Ref. |§ that the critical exponent v of the 
bulk correlation length of the random Ising model does not depend on the dilution, i.e., 
v = 1 at d = 2 both for diluted and not diluted Ising models. The standard (pertubative) 
RG method predicts the change of the universality class by the quenched randomness. 
We think, this is a false prediction. The fact that the standard RG method provides 
incorrect result is not surprising, since it has been demonstrated (in fact, proven) in 
Ref. that this method is not valid at d < 4. 
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